function kap_s = kap_s_pscc(shat,kap0,phi)
%derivative of selling cost per existing customer wrt selling intensity 
kap_s = kap0*shat.^(phi-1);
end

